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ABSTRACT 

We study how the properties of cosmic voids depend on those of the tracer galaxy populations 
in which they are identified. We use a suite of halo occupation distribution (HOD) mocks in a 
simulation, identify voids in these populations using the ZOBOV void finder and measure their 
abundances, sizes, tracer densities, and dark matter content. To separate the effects of bias 
from those of sampling density, we do the same for voids traced by randomly down-sampled 
subsets of the simulation dark matter particles. At the same sampling density, galaxy bias 
reduces the total number of voids by ~ 50% and can dramatically change their size distribu¬ 
tion. The matter content of voids in biased and unbiased tracers also differs. Deducing void 
properties from simulation therefore requires the use of realistic galaxy mocks. We discuss 
how the void observables can be related to their matter content. In particular we consider the 
compensation of the total mass deficit in voids and find that the distinction between over- and 
under-compensated voids is not a function of void size alone, as has previously been sug¬ 
gested. However, we find a simple linear relationship between the average density of tracers 
in the void and the total mass compensation on much larger scales. The existence of this linear 
relationship holds independent of the bias and sampling density of the tracers. This provides 
a universal tool to classify void environments and will be important for the use of voids in 
observational cosmology. 

Key words: cosmology: observations - large-scale structure of Universe - methods: numer¬ 
ical - methods: data analysis 


1 INTRODUCTION 

Several recent studies have advocated the use of large underdensi¬ 
ties in the matter distribution, known as cosmic voids, which can 
be important tools for cosmology in a number of different con¬ 
texts. These include using voids to test modified gravity scenarios 
(Li et al. 2012; Clampitt et al. 2013; Cai et al. 2014; Zivick et al. 
2014), to probe dark energy (Ryden 1995; Lee & Park 2009; Bos 
et al. 2012; Sutter et al. 2015; Pisani et al. 2015), or to detect the ef¬ 
fects of warm dark matter (Tikhonov et al. 2009; Yang et al. 2014). 
These tests make use of observables such as the abundances, size 
distributions, shapes and radial density profiles of voids. This infor¬ 
mation will be available for very large numbers of voids in future 
large-scale survey data such as Euclid (Laureijs et al. 2011). 

Weak lensing by voids (Krause et al. 2013; Melchior et al. 
2014), and their possible gravitational effects on the CMB via 
the integrated Sachs-Wolfe (ISW) effect have also been studied 
(Granett et al. 2008; Cai et al. 2014; Hotchkiss et al. 2015; Planck 
Collaboration et al. 2015). In principle for rare voids such stud¬ 
ies probe the extreme tails of the distribution of fluctuations of the 
gravitational potential so may provide unique insights. Indeed the 
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case of possible ISW imprints of voids is particularly intriguing, 
since theoretical considerations and simulation results show that 
any such detection would be in strong tension with the standard 
ACDM cosmology (Nadathur et al. 2012; Flender et al. 2013; Cai 
et al. 2014; Hotchkiss et al. 2015; Aiola et al. 2015). 

Several catalogues of voids found in galaxy redshift surveys 
have already been published (Pan et al. 2012; Sutter et al. 2012; Na¬ 
dathur & Hotchkiss 2014), and more will be available from future 
data releases. These catalogues can be used for the various stud¬ 
ies described above. However, since the dark matter distribution is 
not directly observable, the comparison of observations with the¬ 
oretical expectation necessary for precision cosmology requires a 
detailed understanding of how biased galaxies actually trace matter 
underdensities. 

A common theoretical model of voids is that based on the ex¬ 
cursion set formalism (Sheth & van de Weygaert 2004), which can 
be modified to account for the bias of galaxy tracers (e.g. Furlanetto 
& Piran 2006). However, this model assumes that voids can be de¬ 
scribed as spherical underdensities surrounded by higher-density 
walls that have reached shell-crossing. Recent results (Falck & 
Neyrinck 2014; Nadathur & Hotchkiss 2015) show that this as¬ 
sumption does not hold for the most commonly used watershed 
void-finders, and more generally may not hold for any practical 
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void-finding technique. Instead voids in simulation and observa¬ 
tion correspond to a much more general set of density minima of 
varying depths. 

This means that the development of cosmological predictions 
for void observables must be guided by and calibrated on simula¬ 
tion results. To properly account for the effects of bias in the galaxy 
distribution, these simulations should be populated with realistic 
mock galaxy populations, for example by using a halo occupation 
distribution (HOD; e.g. Berlind & Weinberg 2002) model. How¬ 
ever, a number of recent studies of voids (e.g. Chan et al. 2014; 
Zivick et al. 2014; Sutter et al. 2015; Pisani et al. 2015) do not fol¬ 
low this approach. Instead these authors randomly down-sample the 
dark matter particle output of simulations to the same mean number 
density as that of galaxies in various surveys, and use these parti¬ 
cles as tracers of voids. This takes into account effects relating to 
the sampling density of tracers. However, the clustering properties 
of sub-sampled dark matter and galaxies are in general expected to 
differ, which should also affect voids. The justification for such a 
procedure has therefore relied on the results of Sutter et al. (2014), 
who suggest that tracer bias has a smaller effect on void properties 
than sampling density and can therefore be ignored. 

This result is however hard to reconcile with the expected de¬ 
pendence of the halo mass function on the background matter den¬ 
sity on larger scales (Sheth & Tormen 2002). Gottlbber et al. (2003) 
find that the numbers of high- and intermediate-mass halos in voids 
is sharply suppressed, and Alonso, Eardley & Peacock (2015) also 
find a strong dependence of the halo mass function on the density 
environment. This means that properties such as the sizes, abun¬ 
dances and density profiles of voids traced by galaxies in such ha¬ 
los should be expected to be very different to those traced by unbi¬ 
ased dark matter. The picture is further complicated by the fact that 
the halo abundance within voids is also expected to be cosmology- 
dependent (Lee 2012; Neyrinck et al. 2014). 

We therefore revisit this question in this paper. We make use 
of a suite of different HOD models to generate biased galaxy tracer 
populations in an A-body simulation, and also use sets of sub¬ 
sampled dark matter particles with the same number densities as 
controls. We use these sets of tracers to identify voids using the 
ZOBOV watershed void-finding algorithm, and compare their prop¬ 
erties. We show that whereas the mean sampling density determines 
the minimum resolvable void size, at the same sampling density 
differences in the tracer bias produce large differences in the void 
size distribution, especially at large void sizes, as well as affecting 
the total number of detected voids and the densities within voids. 
This underscores the importance of using appropriate HOD models 
in simulations when obtaining predictions for void observables. 

We then investigate the distribution of dark matter within 
galaxy voids. The minimum tracer density within each void, n^iini 
varies significantly over the population and we find that it can be 
used to predict the minimum dark matter density within the void. 
The nature of this relationship however depends strongly both on 
the tracer bias and the sampling density. The stacked average den¬ 
sity profile of voids is also of interest, as it could potentially be a 
sensitive test of cosmology (Cai, Padilla & Li 2014, Barreira et al. 
2015). Some studies (Hamaus, Sutter & Wandelt 2014, Nadathur 
et al. 2015) have also suggested the possibility of a universal den¬ 
sity profile, focussing on the variation of the average profile with 
the mean void size. We show that the average profile also depends 
strongly on rimin- 

An important related characteristic of voids is whether the to¬ 
tal mass deficit within them is compensated or not. This is a key el¬ 
ement in determining the environment within voids (Hamaus et al. 


2014). Voids in which the central density deficit is not compen¬ 
sated out to large distances will correspond to gravitational poten¬ 
tial hills, while over-compensated voids are smaller local density 
minima in larger-scale overdense regions, and will therefore corre¬ 
spond to regions of negative gravitational potential. Distinguishing 
between these two cases is therefore crucial for understanding void 
environments, and in particular for the use of voids in the cosmo¬ 
logical tests described above. 

It has previously been suggested (Ceccarelli et al. 2013; Cai 
et al. 2014; Hamaus et al. 2014) that such a distinction can be made 
on the basis of the void size. Instead we find that the transition 
from under- to over-compensation is dependent on both the void 
size and minimum tracer density within the void, and is not well 
predicted by either property alone. However, we are able to iden¬ 
tify a very simple linear relationship between the enclosed mass 
density contrast A at scales much larger than the void and the in¬ 
tegrated average tracer number density within the void, riavg- The 
existence of this linear relationship is found to be remarkably in¬ 
dependent of the bias or sampling density h of the tracer popula¬ 
tion, with Uavg = h marking the transition between A < 0 (under¬ 
compensated) and A > 0 (over-compensated) for all tracer types. In 
contrast, the void size does not provide such a universal diagnostic 
of void compensation. 

In Section 2 below we describe the properties of the simula¬ 
tion, the construction of HOD galaxy samples and the void-finding 
algorithm. In Section 3 we investigate the properties of these voids 
and lay out our results. We discuss their significance and conclude 
in Section 4. Some details of the stacked density profiles around 
voids and comparisons with previous results are presented in Ap¬ 
pendix A. 

2 METHODS 
2.1 Simulation 

We make use of A-body simulation data from the MultiDark sim¬ 
ulation project (Prada et al. 2012), in particular the MDRl data re¬ 
lease.* MDRl simulates a ACDM cosmological model with param¬ 
eters (Dm,f2yv,fifc,/r,nj,cr8) = (0.27,0.73,0.0469,0.7,0.95,0.82) 
using 2048^ dark matter particles in a cubic box of side 1 fi^'Cpc. 
The mass of each dark matter particle in the simulation is 8.7 x 
10^ and it has a force resolution of 7 /i^'kpc. The initial 

conditions for the density evolution are set at redshift z = 65 using 
the Zeldovich approximation. 

Halos are identified in the simulation using the Bound Density 
Maximum (BDM) algorithm first described by Klypin & Holtzman 
(1997), and as updated by Riebe et al. (2013). This algorithm de¬ 
tects as halos spherical regions around density maxima exceeding 
the virialization density threshold in a top-hat model of the growth 
of density fluctuations and removes gravitationally unbound parti¬ 
cles from the halo. The virial radius of the halo, Ryu-, is defined as 
the radius within which the density contrast Avif(z) defined with 
respect to the background density p(z) exceeds the collapse thresh¬ 
old given by the approximation of Bryan & Norman (1998) - at the 
redshift z = 0 snapshot we use, this corresponds to Ayu- = 360. The 
minimum resolved halo mass in MDRl is 1.7 x 10** 

The creation of HOD mocks and other tracer samples for void 
identification is explained below; however, to measure the true den¬ 
sity at void locations we make use of the underlying dark matter 

' Publicly available at www. cosmoslm. org. 
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Table 1. Properties of the mock HOD samples used in this work. Bias values are approx¬ 
imate and intended as a guide only. 


Sample name 

logiffmm 

^logM 

log Mo 

logM; 

a 

n 

{h^Mpc-^) 

b 

Mainl 

12.14 

0.17 

11.62 

13.43 

1.15 

3.18 X 10-3 

1.3 

Main2 

12.78 

0.68 

12.71 

13.76 

1.15 

1.16 X 10-3 

1.4 

LOWZ 

13.24 

0.62 

14.32 

13.20 

0.93 

2.98 X lO-'* 

2.0 


density field from the full resolution simulation output at redshift 
0. This density field is determined on a 1024^ grid using a cloud- 
in-cell interpolation, and is then smoothed using a Gaussian kernel 
of width equal to one grid cell. The resulting density field has a 
sub-Mpc resolution, which is sufficient given the void sizes in our 
catalogues. 


2.2 HOD mocks and tracer samples 


We populate the halo catalogue with mock galaxies using the HOD 
model of Zheng et al. (2007), which assigns galaxies to a dark mat¬ 
ter halo according to a distribution based on the halo mass M. Cen¬ 
tral and satellite galaxies are treated separately. According to this 
model, the mean occupation function of central galaxies is param¬ 
eterized as 


(A^cen(M)) 


1-ferf 


logM — logM„ 


( 1 ) 


and the number of central galaxies in each mass bin follows a 
nearest-integer distribution. On the other hand the number of satel¬ 
lite galaxies follows a Poisson distribution with 


(A^sat(M)) = (iVcen(M)) 


M-Mq 

M' 


( 2 ) 


Central galaxies are placed at the centre of their respective halos, 
while satellite galaxies are distributed through the halo with radial 
distances from the centre drawn from a random distribution based 
on a fiducial NFW mass profile Navarro et al. (1996, 1997). 

The HOD distribution is thus characterized by 5 parameters: 
Mmin, Mq, M'j, ctiogM and a. We choose these parameter to match 
those determined in the literature from fits to the number density 
and clustering properties of observed galaxy samples. 

Two of our mock samples, designated Mainl and Main2, are 
chosen to match volume-limited luminosity threshold samples from 
the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample, 
with r-band absolute magnitude of —20.5 and —21 respec¬ 
tively (Zehavi et al. 2011), with median redshift z ~ 0.1. As these 
authors use a definition of the halo overdensity Ac = 200 with re¬ 
spect to the critical density which differs from the MultiDark def¬ 
inition above, we rescale the halo masses and virial radii from the 
simulation results using the concentration values provided and as¬ 
suming an NFW profile. The large-scale linear bias b, defined by 
8g = bS where Sg is the galaxy density contrast and 8 that of matter, 
is 1.29 for Mainl and 1.40 for Main2. 

Our third sample is chosen to match the properties of the Low 
Redshift (LOWZ) sample of galaxies from the SDSS BOSS data 
releases 10 and 11, with 0.15 < z < 0.43, and is based on the HOD 
parameters obtained by Manera et al. (2015). Although the LOWZ 
sample in reality has a redshift-dependent mean sample density, 
we choose the parameter values for n = 2.98 x 10^^ fi^Mpc^^ and 


uniformly subsample our mocks down to achieve this mean density. 
Manera et al. (2015) use halos within a simulation generated using 
the PTHalos method within second-order Lagrangian perturbation 
theory (2LPT) rather than full A-body simulations; their minimum 
halo mass is also too large to be able to resolve all populated halos, 
so they assign approximately 7% of galaxies to dark matter particle 
positions rather than halos. These differences mean that our mock 
sample is unlikely to precisely match the clustering properties of 
the observed LOWZ galaxies; nevertheless it provides a sufficiently 
realistic realization of a distribution of galaxies with bias b ^2. 

Table 1 summarizes the parameter values and resulting aver¬ 
age number density n and large-scale bias b of our mocks. 

To complement these mock samples and to serve as controls, 
we also use three sets of randomly subsampled DM particles from 
the simulation output as tracers with which to identify voids. These 
samples are referred to as DM Mainl, DM Main2 and DM LOWZ 
and have the same average tracer densities as the Mainl, Main2 and 
LOWZ HOD samples respectively. However, as the random sub¬ 
sampling does not change the fundamental clustering properties of 
the DM field, these three samples are all unbiased {b= 1). Finally, 
we also use three Poisson point samples, consisting of uniformly 
distributed random points with the same mean density as those of 
the mock galaxy tracers. 


2.3 Void finding and merging 

We make use of the ZOBOV watershed void finding algorithm 
(Neyrinck 2008) to identify voids in each of the tracer samples de¬ 
scribed above. ZOBOV uses a Voronoi tessellation field estimator 
(VTFE) technique to reconstruct the tracer density field from a dis¬ 
crete distribution of particles. It then identifies local minima in this 
field and the watershed basins around them. These basins form a 
non-overlapping set of density depressions or voids. 

At this stage, ZOBOV can also merge neighbouring voids to¬ 
gether according to the watershed principle to form a final hierarchy 
of voids and sub-voids. Depending on the choice of input parame¬ 
ters to control void merging, the resulting final void hierarchy and 
distribution of void sizes can vary widely: at one extreme one ob¬ 
tains percolation, with the entire simulation box or Universe being 
included in the largest void; at the other extreme there is no merging 
at all and the original voids are retained. 

Unfortunately, the choice of how to control void merging is 
largely subjective and different options have been applied in the 
literature. Neyrinck (2008) suggests using the ratio of the lowest 
(VTFE-reconstracted) tracer number density along the void wall to 
the minimum VTFE density at its centre, relative to the same ratio 
for spurious voids in Poisson noise, to determine the “most proba¬ 
ble extent”. An alternative is to stop the growth of a void through 
merging if the minimum VTFE density along the watershed ridge 
separating it from a shallower potential sub-void, niink. is greater 
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Table 2. Strategies for controlling merging in the void hierarchy 


Label 

Criteria for 

merging: 


Link density 

Density ratio 

VIDE 

«iiiik < 0.2n 

unconstrained 

Minimal 

niink < « 

r<2 

Isolated 

no merging 

no merging 


than some pre-defined threshold. The VIDE void finding toolkit 
(Sutter et al. 2015), which is based on ZOBOV, sets this threshold at 
0.2 times the mean tracer density fi. However, this value of 0.2 lacks 
any theoretical justification even when the tracers are dark matter 
particles in a simulation, because the recovered voids do not in any 
case correspond to the model objects described by Sheth & van de 
Weygaert (2004) from which the canonical value for the enclosed 
density contrast A = —0.8 derives (Nadathur & Hotchkiss 2015). ^ 

A further difficulty with the VIDE threshold for merging is that 
it applies to the minimum tracer number density along void walls. 
Therefore, when applied to differently biased tracer samples, the 
same numerical value for the threshold corresponds to rather dif¬ 
ferent dark matter densities, and thus to different physical criteria 
for merging. Only for the special case where the threshold value is 
set equal to the mean sample density, i.e. neighbouring voids are 
merged together if < n, does this strategy result in the same 
physical consequences independent of the tracer bias. 

In this paper we will compare three different strategies to con¬ 
trol the merging of voids in the formation of the final hierarchy. The 
first strategy is that applied in the VIDE algorithm, which has been 
used in several recent works (Hamaus et al. 2014; Chan et al. 2014; 
Zivick et al. 2014; Sutter et al. 2015; Pisani et al. 2015; Pollina 
et al. 2015) studying void properties. The second strategy is simply 
to use the original output of density minima from ZOBOV without 
any merging at all, which provides a set of non-overlapping voids 
defined purely by the topology of the tracer density field. We shall 
refer to the voids thus obtained as ‘Isolated’ voids. 

The third strategy we consider allows merging of voids if the 
minimum tracer density on the separating ridge is below the mean. 
However, as this may still allow genuinely distinct shuctures to be 
classified as a single void, we implement an additional condition 
on the density ratio described above, i.e. we prevent the merging of 
two voids if the minimum density on the watershed ridge is more 
than twice the minimum density of the shallower void. This choice 
is based on the Poisson noise characteristics discussed by Neyrinck 
(2008): it is designed to count neighbouring voids as distinct if the 
ridge separating them could not be an artefact of the discrete sam¬ 
pling of the density field. We refer to this as the ‘Minimal’ choice 
of void merging criteria; it has also previously been used by Na¬ 
dathur & Hotchkiss (2014); Hotchkiss et al. (2015); Nadathur et al. 
(2015). The three strategies for void merging are summarized in 
Table 2. 

Most qualitative conclusions regarding the properties of voids 
in the simulation are independent of the particulars of the merging 
strategy chosen; some exceptions to this statement are highlighted 
below. However, quantitative results will inevitably depend on this 


choice. Since the ‘Isolated’ option of preventing merging altogether 
provides a void definition that depends only on the topology of the 
tracer density field without further arbitrary parameters, it may be 
the one most easily described by a future theoretical model. There¬ 
fore, unless otherwise stated, the results presented in this paper re¬ 
late to this case. 


2.4 Defining void properties 

Based on the operational choices described above, we obtain a cat¬ 
alogue of voids within each of our tracer samples. Each void con¬ 
tains a number of tracer particles, and the total void volume is the 
sum of the volumes of the Voronoi cells of each of its member 
particles in the tessellation. We define the void effective radius, 
Rv = (3V/47r)*/^, to be the radius of a sphere with the same vol¬ 
ume as the void, although note that individual voids may be highly 
aspherical. The void member particle associated with the minimum 
VTFE-reconstructed tracer density (i.e., the largest Voronoi cell) 
is labelled the core particle, and this density is denoted by nmin- 
We previously showed (Nadathur & Hotchkiss 2015) that for sub¬ 
sampled dark matter tracers nmin is tightly correlated with (but not 
the same as) the minimum dark matter density within the void, and 
inversely correlated with R^, so that larger voids contain deeper 
density minima. 

The centre of each void is defined as the point of intersection 
of the Voronoi cells of the void core particle and its three lowest- 
density mutually adjacent neighbours in the tessellation, which is 
the point within the void maximally distant from any tracer parti¬ 
cles. Compared to the alternative volume-weighted barycentre def¬ 
inition often used in the literature, this definition is much better at 
locating the true dark matter density minimum within voids (Na¬ 
dathur & Hotchkiss 2015). Unlike the barycentre definition, it is 
also independent of the details of void merging and is less affected 
by shot noise. 

Given the total number of member tracer particles of a void 
and its total volume, it is also simple to define an integrated average 
tracer density for the void, navg- This density is obviously always 
greater than nmini it is also well known (e.g. Nadathur & Hotchkiss 
2014; Nadathur et al. 2015) that the operation of the watershed al¬ 
gorithm means that for some voids rzavg can be much larger than 
Wtnin and even > n due to the inclusion of higher-density filaments, 
walls and halos within the watershed basin constituting a void. 

The location of the centre together with nmirn ”avg and R^ con¬ 
stitute the core observable properties of a void, which can be de¬ 
fined from the tracer population alone. ^ For the same set of minima 
in the tracer density field, navg and Ry will differ depending on the 
details of the void merging criteria described above, but the location 
of the void centre and are independent of this choice. 

Using the gridded dark matter density field from the simu¬ 
lation output we can also measure the true density contrast 8 = 
p/p—l within the void and define the total enclosed density con¬ 
trast within a sphere of radius r, centred at the void centre: 

A(r) = 4 

Jo 


p{y) 


- 1 


y^dy. 


(3) 


^ Even if such a correspondence existed between the theoretical model and 
practical void-finding algorithm, there is no clear relationship between niink 
and A — note that a restriction on nwnk < does not restrict the average 
density of the void to be less than 0.2 — so this choice of threshold value 
would still be unmotivated. 


^ To this list one could also add some measure of the asphericity and direc¬ 
tion of alignment of the axes of aspherical voids, but we will not consider 
these properties in this work. The ratio of tracer densities at the void cen¬ 
tre and walls and the full tracer density profile around the void are also in 
principle measurable, but can be complicated by survey boundary effects in 
observational data (Nadathur & Hotchkiss 2014; Nadathur et al. 2015) 


© 0000 RAS, MNRAS 000, 1-13 







The nature of voids 11 5 



l.Op 

0.8 

0.6 

0.4 

0.2 

0 . 0 - 

0 



DM Mainl 


S= 3 . 18 - 10 “® 


1.0 

0.8 

0.6 

Jo.4 

0.2 


80 


O.OL 

0 



Mainl 


18 - 10 “^ 


80 


R, [/i-^Mpc] 








R, [/f-^Mpc] 


R^ [/i-^Mpc] 

1.0 

DM LOWZ 

1.0 

LOWZ 

1 ' 

1 ' 

0.8 

' .. n = 2.98 10 “^ 

1 

0.8 

. ^= 2.98 10 ^“ 

0.6 

K\ 

In 

O 

K\ "■''■■-■i:;--.. 

0.4 


Jo .4 


0.2 

\ . 

\ 

0.2 

\ “ ”lBf '. 

\ 

\ 

0.0 

_ ^^ _ 

0.0 

_ ^^ _ 


0 20 40 60 80 ^'0 20 40 60 80 


R^ [/i-^Mpc] 





n 

o 

c 

p 


Figure 1. Binned distribution of void minimum tracer densities and radii for all tracer populations as indicated. Plots in the left-hand column are for 
unbiased subsampled dark matter particle tracers, and those in the right-hand column are for the corresponding galaxy tracers with the same average number 
density «. Dotted contours in each plot enclose 95% and 99% of all spurious voids found in uniform random distributions of points with the same n but 
no clustering. The dashed line indicates the minimum resolvable void size as a function of density, and the arrow shows the approximate mean interparticle 
separation . Note that the number of bins remains the same in each plot, so bins cover larger ranges in the right column. 


Within the void, at r < R^, values of A vary widely for voids of dif¬ 
ferent radius and (Nadathur & Hotchkiss 2015). In this paper, 
a property we consider is the compensation of the void underden¬ 
sity on scales much larger than the void, for which we will use the 
value of A within a sphere three times the void radius, A(r = 3Rv), 
as a diagnostic. 

3 PROPERTIES OE VOIDS 

3.1 Effect of tracer bias on void abundances and sizes 

Figure 1 shows the distribution of void sizes R^ and core tracer 
densities for voids each of the tracer samples described in Sec¬ 
tion 2.2. The right hand plot in each row shows the distribution for 
one of the HOD mock samples (Mainl, Main2 or LOWZ); on the 


left is the corresponding distribution for the voids found using the 
dark matter particles subsampled down to the same average tracer 
density n. Therefore within each row, the sampling density remains 
the same but the bias increases. In the left-hand column, the bias re¬ 
mains the same (i.e. b=\) but the sampling density decreases. The 
dotted contours on each plot show the 95% and 99% C.L. contours 
for the equivalent distributions for voids found by the same algo¬ 
rithm in the random Poisson point samples: these contours there¬ 
fore remain the same from left to right in a row. The dashed curve 
in each plot indicates the minimum resolvable void size at the given 
n (Nadathur & Hotchkiss 2015). To clarify the discussion below, for 
the populations shown in this figure no void mergers have been al¬ 
lowed (corresponding to ‘Isolated’ voids described in Section 2.3). 
However, the qualitative features of the figure do not depend on this 
choice. 

Several effects of tracer bias can be discerned from this figure. 
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Figure 2. The cumulative number density of voids with radius > R,,, for all voids in the unbiased DM Mainl sample, and for those in the Mainl HOD galaxy 
mocks. The mean tracer density n is the same for these two samples, so differences in the void sizes are purely due to galaxy bias. The three panels correspond 
to three choices of the criteria to handle void merging described in Section 2.3 and Table 2. The bias of the HOD mocks decreases the total void abundance by 
~ 50% in all cases, and has a large effect on the shape of the distribution, especially for the merging criteria implemented in the VIDE toolkit. 


As discussed in Nadathur & Hotchkiss (2015), the characteristic 
‘banana’ shape of the distribution is naturally common to voids in 
all samples. However, at a quantitative level void sizes and core 
tracer densities depend strongly on the tracer bias as well as the 
mean tracer density. Although the minimum void size is determined 
by the resolution limit set by the sampling density n, voids in mock 
galaxy samples span a larger range of sizes than those in corre¬ 
sponding sub-sampled dark matter populations. Even more impor¬ 
tantly, galaxy voids are emptier of tracers than dark matter voids, 
as evidenced by the distribution of minimum tracer densities nmin 
extending to lower values in the right column. 

This difference in distributions has important conse¬ 
quences for void merging if the criterion for allowing merging is 
based — as in the VIDE case discussed in Section 2.3 — on an ar¬ 
bitrary number density threshold. To see this, note that according to 
the VIDE choice of merging criteria, a necessary but not sufficient 
condition for two neighbouring voids to be merged together as part 
of a larger void is that both must have < 0.2h. As can be seen 
by comparing columns of Figure 1, a much larger fraction of galaxy 
voids satisfy this condition than dark matter voids. This means that, 
if merging is controlled by such a simple threshold alone, void 
mergers will be much more common for voids in galaxy distribu¬ 
tions than in sub-sampled dark matter. Similarly, mergers would 
be more common between Mainl voids than those in LOWZ. This 
will necessarily increase the discrepancy between the observable 
properties of voids in the two populations, as we show later. 

It is also noticeable that at each sampling density n the dis¬ 
tribution of voids in sub-sampled dark matter populations shows 
a greater overlap with that of spurious Poisson voids than that of 
the corresponding voids in biased galaxy tracers. An intuitive ex¬ 
planation for this is that random sub-sampling increases shot noise, 
obscuring the true clustering properties of dark matter and increas¬ 
ing spurious void detections. For galaxy samples, the decrease in n 
is accompanied by increasing bias and therefore enhanced cluster¬ 
ing, so they should be less affected by shot noise. However, we 
caution that just because a void identified in the simulation has 
values compatible with those found for spurious Pois¬ 
son voids it does not necessarily follow that this void is itself an 
artefact of noise! While results for individual voids vary, we find 


that statistically speaking voids within the overlap region still tend 
to correspond to true matter underdensities. For this reason we do 
not apply any of the conservative cuts to the void population previ¬ 
ously used in the literature (see Nadathur & Hotchkiss 2015 for a 
discussion). 

Finally, we also find that the overall abundance of voids de¬ 
creases significantly when the tracer bias is increased at fixed sam¬ 
pling density Ji. In fact for each of the sampling densities consid¬ 
ered, the total number of voids in unbiased sub-sampled dark matter 
particle populations exceeds that in biased galaxy populations by 
~ 50%. This is significantly larger than the variation in void abun¬ 
dance due to changes in cosmological model described by, e.g., 
Zivick et al. (2014); Sutter et al. (2015); Pisani et al. (2015). We 
conclude that predictions for this key observable that are obtained 
from calibration with simulations that do not include realistic mock 
galaxy populations are not observationally relevant. 

To demonstrate this difference in void abundances more 
clearly, in Figure 2 we plot the cumulative number function of 
voids as a function of their size for the DM Mainl and Mainl 
HOD samples, which both have the same sampling density n = 
3.18 X 10^^fi^Mpc^^. Since the three choices of void merging 
strategy produce different results, they are shown separately, as in¬ 
dicated by the titles of each panel. The overall deficit of galaxy 
voids is clearly visible and is due to a lack of small-to-intermediate 
sized voids, in the range 5 <Rv < lOh^^Mpc. 

At large R^,, the number density of galaxy voids is strongly 
dependent on the void merging criteria. The biggest difference be¬ 
tween the DM-only and HOD distributions at large Ry is obtained 
for the VIDE choice of merging threshold, for the reasons explained 
above. Note that the VIDE toolkit has been used by a large number 
of recent void studies which describe the properties of voids using 
down-sampled dark matter tracers as a substitute for galaxies, e.g. 
Sutter et al. (2014); Chan et al. (2014); Zivick et al. (2014); Sutter 
et al. (2015); Pisani et al. (2015). In fact, for the Mainl HOD sam¬ 
ple the VIDE criterion leads to percolation through the simulation 
box, since the largest void has R^, > 500/i^*Mpc. Such percola¬ 
tion only appears when void merging is controlled using a simple 
density threshold alone, as is the case for VIDE. It is strongly de¬ 
pendent on the choice of the arbitrary density threshold, as well as 
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Figure 3. The cumulative number density of voids with radius > S,,, for all voids in the three different HOD galaxy mocks. The panels correspond to the 
different choices of the criteria to handle void merging described in Section 2.3. For the ‘Isolated’ and ‘Minimal’ merging criteria sparser and more highly 
biased galaxies consistently trace larger voids, as expected, but the situation is reversed for the VIDE choice of void merging. 


the sampling density and bias of the tracers. Such percolation may 
generally be considered an undesirable feature in choosing a void 
definition. 

Note that since the topology of the tracer density field is bias- 
dependent, even for the case of ‘Isolated’ voids with no merging the 
void size distribution differs between HOD mocks and dark matter 
tracers at the same sampling density. For the ‘Minimal’ merging 
criteria «!,(> Rf is independent of tracer bias for large but sig¬ 
nificant differences remain in the distribution at small and interme¬ 
diate scales, and in the total void abundance. 

Figure 3 shows the differences in the void size distribution 
over the three mock galaxy samples. Intuitively, increasing tracer 
bias should shift the distribution to larger void sizes as well as de¬ 
creasing the total number of voids: more highly biased galaxies live 
in higher mass halos and are therefore confined to higher-density 
environments, resulting in more large voids, while the lower sam¬ 
pling density n inhibits the ability to resolve small voids. Indeed 
this is what is seen in the middle and right panels of Figure 3, for 
the ‘Isolated’ and ‘Minimal’ merging choices. However, this intu¬ 
itive picture is inverted for the VIDE case, for which the maximum 
void size in each sample instead decreases with increasing galaxy 
bias. This is because of the fixed threshold on tracer number den¬ 
sity for merging, which has different physical consequences in each 
galaxy sample. It provides another reason to disfavour this choice 
of merging. 


3.2 Dark matter within voids 

Although it is the dark matter content of voids that is of primary 
interest for cosmology, in practice it is only the tracer number den¬ 
sities at void locations that can be directly observed. These two 
quantities are obviously distinct for galaxy tracers, but Nadathur 
& Hotchkiss (2015) showed that the tracer number density is not 
the same as the true matter density even when the tracers are 
themselves a random subset of all dark matter particles. This is 
especially true when the number densities are measured by sim¬ 
ply counting tracer particles in a volume around special locations 
within the void, such as the void centre. In fact, since as described 
in Section 2.4 the void centre is by construction far from all trac¬ 
ers, naively reconstructing the tracer number density within a small 


enough volume around the centre returns a tracer number density 
contrast = — 1 for practically all voids. As the hue matter densi¬ 
ties within voids vary quite widely, this means that this naive mea¬ 
sure of the tracer number density is not useful. 

However, the VTFE-reconstructed tracer density is a better in¬ 
dicator of the true matter density (Nadathur & Hotchkiss 2015). 
This is because it is less subject to sparse sampling effects (in¬ 
deed this is the main reason why the VTFE density estimator is 
used in the first place). We find that the minimum or core VTFE- 
reconstructed tracer density contrast within the void, = 

ttmin/h— 1 , is a good predictor of the underlying matter density 
contrast Smin at the void centre.^ Figure 4 shows measurements of 
5min for voids binned by their values of 1 , for 

each of our tracer samples. There is a simple linear relationship be¬ 
tween the two, as evidenced by the goodness of the straight-line fits 
to the data. However, note that 7 ^ Smin even when using dark 

matter tracers, and 7 ^ b8^{^ for voids in the galaxy mocks. 

Some qualitative aspects of Figure 4 can be understood in 
terms of the role of bias and subsampling. At low values of 
voids in galaxy mocks on average correspond to shallower dark 
matter underdensities than do voids in the subsampled dark mat¬ 
ter tracers with the same This is a result of the overall bias 

of the galaxies and the suppression of higher mass halos in low 
density regions — such regions therefore look even emptier in the 
galaxy number density. However, as increases towards zero, 

the situation reverses: this is because subsampling of the dark mat¬ 
ter particles increases the effect of shot noise and means that shal¬ 
low depressions in the tracer density are less likely to correspond 
to true matter underdensities in this case than for the biased galaxy 
mocks. This is also consistent with the picture from Figure 1, where 
voids in subsampled dark matter are much more likely to overlap 
with those in equivalent Poisson point sets. 

These qualitative features are common to voids in each of the 
Main 1, Main2 and LOWZ samples. However, at a more quantitative 

Strictly speaking the correspondence between the void centre and the 
dark matter density minimum is statistical and may not hold for an individ¬ 
ual void. However, any other location within the void that can be deduced on 
the basis of the tracer positions alone will on average always have a higher 
matter density. 
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Figure 4. The binned average dai'k matter density at the void centre as a function of the minimum tracer number density determined from the VTFE recon¬ 
struction, for voids in different tracer populations. In each plot, green circles are for voids traced by sub-sampled dark matter particles and red squares are for 
voids in the corresponding biased galaxy mocks. Error bars indicate the lc7 errors in the mean. The best straight-line fits to the unbinned data are shown with 
dashed (dotted) lines for dark matter (HOD) voids in each case. At low densities, voids in HOD galaxy mocks show smaller values of due to biasing 

effects. Shallow voids in the galaxy distribution are less affected by shot noise due to subsampling and so trace deeper matter underdensities. 


level a cursory search did not reveal a simple way to relate the slope 
and intercept of the straight-line fits to the tracer bias and mean 
number density. This suggests that in order to assess the depth of 
the matter underdensity in voids in observational data, calibration 
using realistic HOD mocks is required on a case-by-case basis. 

A quantity of significant cosmological interest is the distri¬ 
bution of mass within galaxy voids, which can be determined by 
measuring the spherically averaged stacked matter density profile 
about void locations. Previous studies have generally attempted to 
describe the variation of the stacked profile with the mean size of 
the voids included in the stack. Barreira et al. (2015) have proposed 
a functional form with five free parameters. 


P{k) 

P 


1 + 5c 


\1 + {P/S2f ) ’ 


(4) 


where / = r/R^,, to describe the resulting profiles. In Appendix A, 
we show that this function can provide a satisfactory fit to the 
stacked profile data for all void sizes Ry, and examine similar fitting 
forms proposed in other studies (Hamaus et al. 2014; Sutter et al. 
2014; Nadathur et al. 2015). Unfortunately the behaviour of the 
fitted parameters with Ry differs across the different mock galaxy 
samples, meaning that eq. 4 is not predictive, i.e., the best-fit pa¬ 
rameters must be determined on a case-by-case basis for voids in 
each sample. 

In any case, describing the best-fit profile parameters as a 
function of Ry alone is not a useful strategy, since it obscures the 
wide variation in tracer densities within voids of the same size 
that can be seen in Figure 1. Indeed Nadathur & Hotchkiss (2015) 
showed explicitly that the mass density profile of voids is also 
strongly dependent on the tracer values even for fixed Ry. To 
demonstrate the full variety of watershed void properties, in the 
next section we consider the dependence on njj,in and Ry simulta¬ 
neously. 


3.3 Compensation of mass around voids 

Describing the full variation of the density distribution around void 
centres in terms of a fitting formula such as eq. 4 and the depen¬ 


dence of its free parameters on void properties is a complicated 
task. Instead we choose to describe the density environment around 
the void location by a single number, namely the total enclosed 
mass density contrast A within a sphere of radius three times the 
effective void radius Ry. This radius is large enough that A under¬ 
goes no further zero-crossings at larger r. Therefore the sign of 
A(r = 3Ry) distinguishes between two important classes of void 
environments. Overcompensated voids, with A > 0, are local den¬ 
sity minima embedded within larger scale contracting overdensi¬ 
ties, reminiscent of the void-in-cloud scenario described by Sheth 
& van de Weygaert (2004). On the other hand, for A < 0 the inte¬ 
rior mass deficit of the void is not compensated by surrounding high 
density regions. Such voids will correspond to regions of positive 
gravitational potential, <I> > 0. Being able to distinguish between 
these two cases on the basis of primary void observables is im¬ 
portant for studies of the gravitational effects of voids, particularly 
measurements of lensing and secondary CMB anisotropies due to 
voids. 

Figure 5 shows the dependence of A(r = 3Ry), measured from 
the full resolution simulation output, on and Ry for voids in 
each of the three mock galaxy samples. The same trends are seen 
across all samples, suggesting the possibility of a universal descrip¬ 
tion of the void environment. However, contrary to previous sug¬ 
gestions (Ceccarelli et al. 2013; Cai et al. 2014; Hamaus et al. 2014) 
the transition from over-compensation to under-compensation is 
not a function of the void size alone: both the smallest and largest 
values of A are found at small void radii. Shallower voids are on 
average more likely to be over-compensated, so a cut based on 
the minimum tracer density at roughly ^rnin/” 0-3 (s® applied by 
Hotchkiss et al. 2015) would better distinguish under-compensated 
voids, but is still not completely satisfactory. 

The dependence of void properties on variables other than the 
void radius is significant for their theoretical modelling. Thus far, 
such models have followed Sheth & van de Weygaert (2004), who 
develop an excursion set model roughly analogous to that for ha¬ 
los. In this scenario, the void radius Ry is analogous to the halo 
mass and so all void properties, an in particular, the variation of 
the average density profile, are described purely with reference to 
Ry. Figure 5 shows that this analogy does not work well, because 
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Figure 5. The variation of void compensation as a function of minimum tracer density and void size. The colourbar indicates the binned average values of the 
total matter density contrast enclosed within a sphere of radius 3Rv centred at the void locations. The same trends are seen for voids in all three galaxy mocks: 
small and shallow voids are on average over-compensated, while deep and large voids are under-compensated, but the variation depends on a combination of 
^min and is not satisfactorily accounted for by either variable alone. 
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Figure 6. The mean enclosed matter density contrast A(r = 3i?v) for voids as a function of the void radius R^, alone. Left: Comparison between voids in the 
Mainl HOD sample and those in the unbiased DM Mainl tracers, which have the same mean tracer density n. Similar results hold for other pairs of samples 
with the same tracer density. Right: Voids in the three HOD mock galaxy samples. In all cases, small voids are on average overcompensated (A > 0) and 
large voids on average exactly compensated (A ~ 0). The radius scale for the approach to compensation depends on both h and the tracer bias. Information on 
undercompensated voids seen in Figure 5 is lost when through averaging over «min. 


it fails to account for the strong systematic dependence on Umin- 
In a companion paper (Nadathur & Hotchkiss 2015), we have also 
shown that watershed voids fail to satisfy several other important 
assumptions of the Sheth & van de Weygaert (2004) model. 

It is however possible to ignore the dependence on timin and 
describe the average behaviour of A as a function of alone, as 
has been done in previous works. The resulting behaviour is shown 
in Figure 6. For all samples, small voids are on average overcom¬ 
pensated (A > 0), but approximate compensation is reached at large 
Ry. Note that all information of the undercompensated voids clearly 
visible in Figure 5 is lost in marginalizing over Kmin- The radius 
scale above which voids are on average compensated is dependent 
on both the tracer density and bias. 

In Figure 7 we show the dependence of A instead on the aver¬ 
age tracer number density within a void, Wavg, for all six tracer pop¬ 
ulations used in this study. Using Wavg as the independent variable 
preserves information on both under- and overcompensated voids. 


In all cases a simple linear relationship exists between navg and 
A. Straight-line fits to the data show that Uavg = « is the universal 
dividing line between under-compensated and over-compensated 
voids in every case. We conclude that, in contrast to R-^, ns,vg is a 
universal indicator of the compensation of the mass deficit within a 
void, independent of the bias properties or sampling density of the 
tracer population in which the voids are identified. 

At one level, this result appears intuitively simple: regions 
containing a higher than average number density of tracers, whether 
dark matter particles or galaxies, also correspond to a higher than 
average matter density. However, note that fiavg is defined for voids 
of arbitrary — and generally aspherical — shape determined by the 
topology of the tracer density field, whereas A refers to the mass 
density within idealised spheres. There is also a large difference in 
scales, since A is measured within a volume 27 times larger than 
that of the void. 

It is also worth noting that the relationship between local 
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Figure 7. Void compensation as measured by the enclosed matter density contrast A{r = 3Sv) from eq. 3, as a function of the overall density of tracers within 
the void, Mavg- The left panel shows the trends for voids in suhsampled dark matter, and the right panel for voids in galaxy tracers. The dashed lines show the 
best fits to the data; in all cases a simple linear relationship provides an excellent description. The value Mavg = « provides a universal distinction between 
under-compensated and over-compensated voids, irrespective of the properties of the tracer population in which voids are identified. 


tracer and matter densities at the void centres is not independent 
of the properties of the tracer population, as shown in Figure 4. 
Voids in different samples trace minima of very different depths. 
Similarly we were also unable to find any comparable relationship 
between the locally-measured tracer and matter densities at void 
locations, i.e. the profiles of 8{r) are not proportional to those of 
8n{r) for any tracer type. This simple linear relationship exists only 
for integrated quantities Wavg and A. 


4 CONCLUSIONS 

The existing theoretical model of voids based on the excursion 
set formalism of Sheth & van de Weygaert (2004) has previ¬ 
ously been shown to be inadequate for describing the properties 
of voids that can be practically identified in simulations by means 
of the commonly-used watershed algorithm (Nadathur & Hotchkiss 
2015). This means that meaningful predictions for void properties 
and the cosmological constraints that may be obtained from the use 
of voids can at present only be obtained through a simulation-led 
approach, but to be observationally relevant such simulations must 
aim to replicate the types of voids that would be seen in actual 
galaxy survey data. 

Our aim in this paper has been to provide a comprehensive 
examination of the twin roles of galaxy bias and sparsity in de¬ 
termining the properties of voids, and in the process to improve 
existing treatments in the literature. To do so we made use of a 
suite of mock galaxy samples in an V-body simulation, which were 
used as tracer populations for void identification. These void cat¬ 
alogues were compared with equivalent ones obtained using unbi¬ 
ased subsets of the dark matter particles in the simulation as tracers. 
Although our results were obtained specifically using the popular 
ZOBOV algorithm, the general qualitative features of voids we de¬ 
scribe will be applicable to any void-finder based on the watershed 
principle. 


Our first practical result has been to show that, contrary to 
previous claims (Sutter et al. 2014; Zivick et al. 2014; Sutter et al. 
2015), tracer bias plays a crucially important role in determining 
the overall abundance of voids, their size distributions, and mat¬ 
ter content. The differences in such void properties introduced by 
tracer bias are comparable to or larger than those reported from 
differences in cosmology (Zivick et al. 2014; Yang et al. 2014; Sut¬ 
ter et al. 2015; Pisani et al. 2015). This result is not entirely unex¬ 
pected, since all of these quantities naturally depend strongly on the 
clustering properties of the tracer population, which are in turn re¬ 
lated to the tracer bias. However, it means that the commonly-used 
strategy of approximating galaxy voids by simply down-sampling 
the density of dark matter particles in a simulation is incorrect and 
that predictions for void observables calibrated on such simulation 
results will unfortunately bear little or no relation to data obtained 
from real galaxy surveys. A similar conclusion has been indepen¬ 
dently reached by Pollina et al. (2015). 

Our results have also clarified the general nature of watershed 
voids. Such voids correspond to density minima of widely varying 
depths typically in the linear regime, rather than the highly non¬ 
linear evolved objects postulated by some theoretical models. Most 
previous studies of voids, from both theoretical and simulation per¬ 
spectives, have focussed on describing void properties as functions 
of the void size alone. In reality the properties of voids of the same 
size show a strong systematic dependence on other factors as well. 
Our work here and in the companion paper Nadathur & Hotchkiss 
(2015) has aimed to highlight the important role of tracer number 
densities in particular as predictors of the matter content of voids. 

From this perspective, the most important results in this paper 
are shown in Figures 5 and 7, which show how the large-scale mat¬ 
ter fluctuation around void locations can be predicted from easily 
measurable properties of the tracer galaxies within the void, irre¬ 
spective of the detailed clustering properties of the galaxy popu¬ 
lation. In particular, the average galaxy number density within the 
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void can serve as a universal diagnostic of the mass compensation 
on much larger scales. 

The nature of this mass compensation is interesting in its own 
right and as a clue towards the formation and evolution history of 
voids. However, it is primarily important for the use of voids in 
cosmology because it serves as a signifier of the value and sign of 
the gravitational potential <I> at the void location. For instance, for 
any observational measurement of the gravitational effects of voids, 
such as through lensing distortions or secondary CMB anisotropies, 
one would at a minimum like to exclude voids corresponding with 
•F < 0. In the past, various authors have attempted to use other di¬ 
rectly observable properties of watershed voids, such as the ratio of 
densities in the void wall and at its centre (Granett et al. 2008), the 
void size (Cai et al. 2014; Planck Collaboration et al. 2014), or the 
minimum galaxy number density within the void (Hotchkiss et al. 
2015) as proxies for the value of <F. Our result indicate instead that 
navg is a better indicator, and that this holds for voids seen in any 
galaxy sample. A generalization of this measure to a weighted av¬ 
erage density for galaxies with a redshift-varying selection function 
should also be quite straightforward. 

We should at this point introduce a note of caution. Strictly 
speaking, our results in this paper concern the compensation of 
mass in voids, and not the potential itself. We have used A as a 
single number with which to characterize the total large-scale mass 
deficit at void locations, but the gravitational potential depends also 
on the details of the density profile and the size of the void, as it is 
related to an integral over r of A(r). Nevertheless, it remains the 
case that n^yg is an excellent predictor of <F, as we will show in 
forthcoming work. 
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Figure Al. Average density profiles for voids of different sizes. Voids in each mock galaxy population are binned on the basis of their radius R„ and stacked 
so that their centres coincide. The average matter density distribution for the stack is then determined from the dark matter field of the simulation. Solid lines 
show the best-fit forms of the profile in eq. 4 in each case, and dashed lines the best-fit forms of eq. Al. Some radius bins are omitted for clarity. 


APPENDIX A: VOID DENSITY PROEILES AND FITTING 
FORMULAE 


AI Matter density profiles 

Although stacked void density profiles have been studied by several 
authors, results in the literature are rather contradictory (e.g Riccia- 
rdelli et al. 2014; Hamaus et al. 2014; Sutter et al. 2014; Nadathur 
et al. 2015; Nadathur & Hotchkiss 2015; Barreira et al. 2015). We 
start by considering the dependence of the average density profile 
on the void radius Rv To do this, we create stacks of voids in bins 
of Rv and centred on the void centre, and determine the average 
profile p{r)/p for each stack using the gridded density information 
from the full resolution simulation output. From this data we deter¬ 
mine the best-fit parameters 8c, S 2 , a and j3 in eq. 4 for each 

stack. 

Figure Al shows the profile data for voids of different sizes 
in each galaxy sample, along with the best-fit forms of eq. 4 (solid 
curves) in each case. Hamaus et al. (2014) have proposed a ‘univer¬ 
sal’ fitting form for the profile, which corresponds to a special case 
of eq. 4 with i2 = 1: 


P \l + {r/R,.fJ 


(Al) 


The best-fit versions of this constrained profile are shown with 
dashed lines in Figure Al. It can be seen that they generally pro¬ 
vide somewhat worse fits to the data, both inside the void and in the 
compensating shell. 

A key feature of both fits is that, as shown in Figure A2, within 
the sample of voids from each population the fitted parameter dc is 
a decreasing function of R,,. This is consistent with (and can be in¬ 
directly deduced from) the generic watershed behaviour shown in 
Figures 1 and 4, and also agrees with the results of Nadathur & 
Hotchkiss (2015) where we found that larger voids contain deeper 
density minima. However, it is opposite to the behaviour claimed 
in Hamaus et al. (2014), who find Sc increases with R,,- For the 
other fitted parameters we find no consistent trends across the dif¬ 
ferent galaxy samples, so we are also unable to confirm the claim 
of universality of the void density profile made by those authors. 
In particular, the dependence of a and /3 on scale radius does 
not follow the pattern described by Hamaus et al. (2014) and Sutter 
et al. (2014) (see Nadathur & Hotchkiss 2015 for more details). 


It should be noted that Hamaus et al. (2014) proposed eq. Al 
as a description of the average density profile for stacks centred 
about void barycentres, whereas in this paper we define the void 
centre differently to better identify the true location of the under¬ 
density. In Nadathur & Hotchkiss (2015) we showed that for some 
high-density tracer populations and for some choices of merging 
criteria, the average density at the barycentre can show a non¬ 
monotonic behaviour with increasing void size (see Figure 3 in par¬ 
ticular). This occurs due to the displacement of the barycentre lo¬ 
cation from the true underdensity in voids composed of very large 
numbers of merged sub-voids. Even in these situations, however. 
Sc does not increase monotonically with R,, as claimed by Hamaus 
et al. (2014). This unusual behaviour also does not apply for any 
of the galaxy samples used in this work, which all show monotoni¬ 
cally decreasing Sc{Rv) irrespective of the choice of void centre. 

However, it is worth stressing again that analyzing the varia¬ 
tion of the mean density profile as a function of R,, obscures the 
important sytematic variation across the void population with the 
values of shown in Figure 5 and discussed in Nadathur & 
Hotchkiss (2015). The stacked profiles in R,, bins in this Appendix 
are therefore included only to enable comparisons with previous 
results in the literature. For voids of the same size, the variation of 
the stacked density profile in stacks of the same is large, as ex¬ 
pected from Figure 5. If voids are stacked according to their or 
Wavg values rather than R,,, the form of eq. Al is completely unable 
to fit the resulting profiles. However, the form of eq. 4 continues to 
provide a good description in all cases. 


A2 Galaxy number density profiles 

In a previous work (Nadathur et al. 2015), we found that the stacked 
galaxy number density profile for a selected subset of voids in sev¬ 
eral different galaxy populations drawn from data from SDSS sur¬ 
veys at different redshifts were very similar, and argued that this 
might be a sign of universality. However, the size of the data sam¬ 
ples used in this work were small, and the errors in determination 
of the profiles were correspondingly large. The number of voids in 
our simulation samples exceeds that used in the previous work by a 
factor of ~ 10: this is because tbe simulation volume is much larger 
than the survey volume for some of the observational samples, and 
boundary effects due to survey masks are not present. We are there- 
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the surviving voids in each of the Mainl, Main2 and LOWZ mock 
galaxy samples separately, without binning by radius, to determine 
the average galaxy number density profiles. For consistency, we 
also use the ’VTFE’ estimator for measuring galaxy density pro¬ 
files discussed in that work and centre stacks on the void barycen- 
tres, though the results are not strongly dependent on this choice. 
The resulting profiles are shown in Figure A3. The results in the 
different galaxy populations are similar to each other, but distin¬ 
guishable with high statistical significance due to the small error 
bars. However they are also consistent with the profile seen in Na- 
dathur et al. (2015) within the much larger statistical errors of that 
study. 


Figure A2. The dependence of the best-fit values of the parameter 5^ in efi- 4 
as a function of the mean void radius of the stacks shown in Figure Al. in 
each case 5^ decreases with and this trend is unchanged if the alternative 
fitting form in eq. Al is used, or for stacks centred on the void barycentres. 
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Figure A3. Galaxy number density profiles for voids in the three mock 
galaxy samples. To match the methodology used by Nadathur et al. (2015) a 
selection cut «niiii/H < 0-3 is applied and stacks are centred on void barycen¬ 
tres. The resulting average profiles in the different galaxy populations are 
similar but statistically distinguishable given the large sample size. The dot¬ 
ted lines are linear interpolations and do not represent fits to the data. 


fore able to re-examine the universality of galaxy density profiles 
with greater statistical precision. 

Nadathur et al. (2015) analysed only those voids with 
Timin/^ < 0.3. In the absence of information of the true matter con¬ 
tent of galaxy voids, this very conservative cut was intended to 
remove Poisson contamination of the void sample, based on cal¬ 
ibration of the probability P(«tj,i„|Poisson). Ad hoc cuts to the void 
population on the basis of void radius Ri, have been advocated else¬ 
where in order to achieve the same purpose. Our analysis of the 
matter density content of voids in simulation, both in this paper 
and in Nadathur & Hotchkiss (2015), indicates that on average all 
voids correspond to true matter underdensities, so such cuts are not 
necessary and lead to sub-optimal use of available data. However, 
for the purposes of direct comparison with the results of Nadathur 
et al. (2015), we apply the same methodology here. We then stack 
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